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ABSTRACT 

3C 318 is a z — 1.574 radio-loud quasar. The small physical size of its radio jets indicate that these jets 
were triggered relatively recently. In addition to the ultraviolet continuum emission being reddened by 
dust, detections with IRAS and SCUBA show it to have an exceptionally high far-infrared luminosity. 
We present CO (2-1) observations of 3C318 made with the IRAM Plateau de Bure Interferometer. 
We detect CO (2-1) emission with a FWHM = 200 km s" 1 at a signal-to-noisc ratio of 5.4. There is 
evidence for positional (~ 20kpc) and velocity (~ — 400kms~ ) offsets between the molecular gas and 
the quasar which may be due to the quasar experiencing a major merger. The mass of molecular gas 
inferred from our observations is Mh 2 = (3.0 ±0.6) x 10 10 Mq. This molecular gas mass is comparable 
to that in sub-mm-selected galaxies at similar redshifts. The large molecular gas mass is consistent 
with the primary source of heating for the cool dust in this quasar to be massive star formation with 
a star formation rate of 1700 M Q yr -1 and a gas depletion timescale of 20Myr. Our observations 
support the idea that star formation episodes and jet triggering can be synchronised. 
Subject headings: quasars individual (3C318) — radio lines: galaxies — galaxies:starburst 



1. INTRODUCTION 

Quasars are believed to be actively accreting black 
holes hosted by massive galaxies. The discovery of su- 
permassive black holes in the nuclei of all nearby galaxies 
with stellar bulges (Magorrian et al. 1998) indicates that 
such activity is not exceptional but the norm. What is 
still unknown is exactly what physical processes are re- 
sponsible for the black hole mass - stellar bulge mass 
correlation. 

Of particular interest is to observe the simultaneous 
growth of black holes and stellar mass in galaxies at the 
peak epoch of activity (redshifts 1 < z < 3). However, 
surveys of sources detected in hard X-rays (probing AGN 
activity) and sub-millimetre emission (probing star for- 
mation) found very little overlap amongst the brightest 
sources in each class (Almaini et al. 2003; Waskett et al 
2003). Only a low level of black hole accretion is found for 
most sub-millimetre-selected galaxies (SMGs; Alexander 
et al. 2005). Therefore the most intense periods of star 
formation and black hole growth are typically decoupled. 
However, to be consistent with the black hole mass - stel- 
lar bulge mass correlation, models have been proposed 
whereby the X-ray bright and sub-mm bright phases are 
part of an evolutionary sequence within a single event, 
such as a major merger (Granato et al. 2001; Archibald 
et al. 2002). 

Radio-loud quasars possess a particularly useful fea- 
ture. Their relativistic jets can be modelled to infer the 
age of the radio source and hence the time when the 
jet was triggered. Recent multi-wavelength observations 
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have revealed several interesting properties of young ra- 
dio sources: (i) strong associated C IV absorption (Baker 
et al. 2002); (ii) high sub-mm luminosities (Willott et al. 
2002); (iii) fast, massive outflows of neutral gas (Mor- 
ganti et al. 2005). These observations point towards a 
phase of great activity in the host galaxies of sources with 
recently triggered radio jets. 

A long standing issue in relation to the sub-mm detec- 
tions of luminous quasars is whether the radiating dust 
is heated by young stars or the active nucleus (Almaini 
et al. 1999). The correlation between quasar optical lu- 
minosity and sub-mm luminosity (Willott et al. 2003) 
could be caused by direct heating of the dust by the 
active nucleus, or it could instead be due to correlations 
between black hole and galaxy mass (Omont et al. 2003). 
A significant number of sub-mm-luminous quasars have 
been found to possess large reservoirs of molecular gas 
via CO line observations (see Solomon & Vanden Bout 
2005 for a compilation). Although time consuming work, 
such observations are critical because molecular gas pro- 
vides the fuel for star formation. The existence of a large 
reservoir of molecular gas provide useful insights on the 
evolutionary state of the system via its gas-to-dust ratio. 

There are now 16 z > 1 quasars with measured CO 
emission of which at least 11 are known to have their 
CO detections assisted by gravitational lensing. How- 
ever, only one of these quasars, Q0957+561, lies at red- 
shifts 1 < z < 2 when the bulk of the accretion in optical 
quasars occurred (Barger et al. 2005). This redshift 
range contains very few known CO emitters with the 
only other ones being the extremely red starburst galaxy 
HR10 (Andreani et al. 2000) and a gravitationally-lensed 
SMG SMM J02396 (Greve et al. 2005). It is also worth- 
while to note that most CO line observations of high- 
redshift objects are of the high-J transitions (Solomon 
& Vanden Bout 2005). The low- J transitions are less bi- 
ased by the details of the excitation and are more likely 
to offer a reliable estimate of the total molecular gas mass 
(Papadopoulos & Ivison 2002; Hainline et al. 2006). 
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We have used the IRAM Plateau de Bure Interfer- 
ometer (PdBI) to search for low-J CO emission in the 
z = 1.574 radio-loud quasar 3C318. As discussed in 
Willott et al. (2000), 3C318 is interesting because (i) it 
is one of the most distant far-IR sources detected with 
IRAS and therefore has a very high far-IR luminosity, (ii) 
it has a high sub- mm luminosity measured with SCUBA, 
(iii) the quasar ultraviolet continuum is reddened by dust 
and (iv) it is a very compact (i.e. young) radio source. 
Here we report the discovery of CO (2-1) emission in this 
quasar and discuss the implications. Cosmological pa- 
rameters of Hq = 71 km s _1 Mpc -1 , f^M = 0.26 and 
Oa = 0.74 (Spergel et al. 2006) are assumed through- 
out. 

2. OBSERVATIONS & RESULTS 

3C 318 was observed in the 3 mm band with the PdBI, 
in the 6D (6 antennae, compact) configuration, on the 
nights of 2002 April 19 and 22. The weather conditions 
were generally excellent, particularly the seeing (0.74" 
and 1.01" on April 19 and 22 respectively) leading to very 
good phase stability. The central frequency was tuned to 
89.60 GHz, close to 89.5641 GHz which is the frequency 
of the CO (2-1) transition at a redshift of z = 1.574 (the 
rcdshift of the narrow [O in] emission line detected in 
near-IR spectroscopy by Willott et al. 2000). The max- 
imum bandwidth of 560 MHz was used, equivalent to 
1875 km s" 1 . 

The observations lasted a total of 13.4 hours, which in- 
cluded observations of the calibrators 3C 273 (beginning 
of the night) and MWC 349 (end of the night) as well 
as regular visits to the phase calibrator 1413+135. The 
data were reduced at IRAM in Grenoble. After ampli- 
tude and phase calibrations, 3C318 was detected very 
strongly at 3 mm. Due to incomplete coverage of the 
uv plane, the beam includes two artificial sidelobes, al- 
though this is not important for a detection experiment 
and 3C318 appears unresolved anyway. 

The continuum flux density of 28.65±0.28mJy is con- 
sistent with a point source convolved with the beam (of 
size 8.05" x 4.32"), with the location of the peak just 0.2" 
from the position of the compact radio source. 3C 318 has 
been detected at 1.25 mm with MAMBO at the IRAM 
30-m Telescope and the relative fluxes of the 1.25 mm 
and 850 /im emission indicate that the 1.25 mm emission 
is dominated by non-thermal synchrotron (Haas et al. 
2006). Therefore, the 3 mm continuum must be entirely 
due to synchrotron from the compact radio source. 

A search of the data-cube using a routine kindly pro- 
vided by Roberto Neri revealed a fairly broad emission 
line close to the location of the continuum emission. This 
was the only significant (5cr) broad peak to be found in 
the entire field-of-view, giving confidence that it is not 
spurious. The continuum flux density contribution was 
then subtracted from the spectrum (in the uv plane), as- 
suming a flat spectrum through the band. A frequency 
window was previously specified around the line, how- 
ever, so the flux of the continuum to be subtracted did 
not include the line. After this continuum subtraction, 
a spectrum was extracted from the data-cube, using a 
polygon matching the shape of the synthesised beam. 

The spectrum was rebinned to 20 MHz (33kms _1 ) 
channels and is shown in Fig. 1. Positive emission is 
visible in nine adjacent channels, corresponding to a full 
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Fig. 1.— Spectrum of CO (2-1) emission in 3C318. The data 
have been rebinned into 20 MHz (33kms — 1 ) channels. The curve 
shows the result of a single Gaussian fit to the line with a FWHM 
of200kms~ 1 . The upper axis shows the redshift corresponding to 
the CO (2-1) transition. 



width zero intensity (FWZI) of 300 km s^ 1 . Integrating 
the emission over these channels results in a line flux of 
1.19 ± 0.22 Jykms -1 . Also shown in Figgis a Gaus- 
sian fit to the emission line. The best fit Gaussian has 
a centre at -403 km s -1 and a FWHM of 200kms _1 . 
The emission is somewhat more boxy than a Gaussian, 
possibly indicating that a multiple Gaussian fit may de- 
scribe the kinematics better, as has been found for sev- 
eral SMGs (Greve et al. 2005). However, our data only 
have a signal-to-noise ratio for the line of 5.4, so a multi- 
component fit would not be well constrained and is not 
attempted. 

A map of the line emission was made by extracting 
the velocity channels containing the line. This map is 
contoured in Fig.[21 The peak of the line emission is 
2.4" west and 0.5" north of the continuum peak. We 
now consider the significance of this offset, given the rel- 
atively low signal-to-noise ratio (SNR) of the emission 
line. The beam shape in the map is an ellipse with Gaus- 
sian FWHM of 8.05" and 4.32" at a position angle of 79° 
east of north. In the one-dimensional case, the la po- 
sitional uncertainty is given by a x = (FWHM/2)/SNR 
(e.g. Downes et al. 1999). For the two axes of the el- 
liptical beam we find a x = 0.75" and <j y = 0.40". Along 
the angle corresponding to the offset of 2.4" west and 
0.5" north, the size of the lcr positional uncertainty is 
0.65" and hence the significance of the positional offset 
is 4cr. We note that given the low SNR and resolution of 
the line, we cannot rule out a considerable fraction of the 
molecular gas existing at the same site as the synchrotron 
continuum. 

In addition to this spatial offset, the velocity peak of 
the molecular line emission is blueshifted by 400 km s -1 
compared to the redshift determined from the narrow 
[O in] emission. The narrow emission lines of a quasar 
are usually a good indicator of the galaxy systemic ve- 
locity (Vrtilek & Carleton 1985). Greve et al. (2005) 
consider the velocity offsets between molecular and ultra- 
violet/optical emission lines in quasars and SMGs. They 
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Fig. 2. — The contours show the CO (2-1) emission from 
3C318 (integrated over channels corresponding to — 550kms — 1 to 
— 270 kms - 1 ). The location of the 3mm continuum peak (domi- 
nated by synchrotron emission from the compact radio source) is 
indicated with a cross. The offset between these positions is signif- 
icant at the 4<r level. The beam size is 8.05" X 4.32" and is shown 
in the box in the lower-left corner. Contours are drawn at 0.7, 1.05, 
1.40, 1.75 mjybeam -1 . 



find that offsets of several hundred kms are common 
for the molecular and Lya lines for quasars and SMGs. 
However, they find very small offsets (< 100 km s -1 ) for 
the molecular and Ha lines in four SMGs. The explana- 
tion for this is that in SMGs, Ha traces the same star- 
forming regions that produce the molecular emission and 
is not so affected by dust obscuration and absorption as 
Ly a. For the case of quasars it is well-known that there 
are large offsets between the broad and narrow emission 
lines (e.g. Richards et al. 2002) and that it is difficult to 
locate the true Lya peak at high redshift. 

How should we interpret the projected spatial (~ 
20kpc) and velocity (~ 400kms~ ) offsets between the 
radio core (situated close to the supermassive black hole 
at the bottom of the galaxy potential well) and the bulk 
of the molecular gas? We note that similar sized off- 
sets have been found for SMGs between the peak of the 
molecular gas and the peak of the if-band stellar contin- 
uum (Tacconi et al. 2006). The most likely explanation 
is that 3C318 is in the early stages of a major merger, 
when such offsets between the two galaxies are common 
(Patton et al. 2000). In this case, we find that the bulk of 
the molecular gas (and star formation) is in the compan- 
ion galaxy to the active radio source. Higher resolution 
observations of the molecular gas and thermal dust emis- 
sion are clearly needed to fully understand the situation 
in this complex system. 

We calculate the CO line luminosity and molec- 
ular gas mass assuming a conversion factor a = 

as in Solomon & Vandcn Bout 
(2005). The line luminosity is L' co = (3.8 ± 0.7) x 
10 10 Kkms _1 pc 2 and the mass of molecular hydrogen 
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Fig. 3. — Plot of Lfir/Lq (oc SFR timescale) vs Lfir for 
ULIRGS (open circles; Solomon et al. 1997), SMGs (filled circles; 
Solomon & Vanden Bout 2005), quasars and radio galaxies (plus 
signs; Solomon & Vanden Bout 2005) and 3C318 (asterisk). The 
right axis shows the molecular gas depletion timescale assuming a 
star formation rate derived from the far-infrared luminosity. 



estimated to be M Ha = (3.0 ± 0.6) x 10 10 Af . These 
values are comparable with the mean of a sample of CO- 
detected SMGs (Greve et al. 2005). We note that there 
is no evidence for gravitational lensing of 3C 3f 8 (Willott 
et al. 2000). 

To compare the mass of molecular gas with the rate 
at which it is being turned into stars we use the far- 
infrared luminosity to estimate the star formation rate. 
Willott et al. (2000) showed that 3C318 has an ex- 
tremely high far-infrared luminosity with significant con- 
tributions from dust warmer than that which is typi- 
cal of SMGs. Since this warm dust component is likely 
heated by the quasar rather than young stars, we calcu- 
late the far-infrared luminosity using only the SCUBA 
850/im flux-density and assuming a cool (T — 41K) dust 
spectrum as in Willott et al. (2002), yielding the relation 
Lfir w 2 x 10 12 (S^o / mJy) L Q . The measured 850/Ltm 
flux-density of 3C318 is 7.78 ± 1.00 mJy (Willott et al. 
2002). Haas et al. (2006) detect 3C318 with MAMBO 
at 1.25 mm with a flux-density of 5.61 ± 1.42 mJy. This 
suggests some fraction of the 850 /im flux is due to syn- 
chrotron and from consideration of the synchrotron spec- 
tral slope we estimate the thermal contribution at 850 ^m 
to be k 5mJy. Therefore the cool dust far-infrared lu- 
minosity is log 10 Lfir = 13.0 Lq. Integrating over the 
stellar initial mass function, this corresponds to a star 
formation rate of 1700 Af Q yr -1 (Kcnnicutt 1998) and a 
corresponding timescale for molecular gas depletion of 
~ 20Myr. By comparison, the compact radio jets ex- 
tend only 7kpc and, assuming a hotspot advance speed 
of < 0.1c, were triggered < 1 Myr before the time of ob- 
servation. 

In Fig.|3 we show how these values for the far- 
infrared luminosity, ratio of far-infrared to CO luminos- 
ity and molecular gas depletion timescale compare with 
all known z > 1 molecular line emitters from the liter- 
ature (Solomon & Vanden Bout 2005) and low-redshift 
ULIRGs (Solomon et al. 1997). This plot shows that 
high-redshift quasars and radio galaxies occupy a similar 
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location to SMGs and that 3C 318 is typical of both of 
these populations. 

The discovery of a large molecular gas mass in 3C 318 
has implications for the power source of the exceptionally 
luminous dust emission. In the case of dust heated di- 
rectly by the active nucleus, a large molecular gas mass 
is not required. The measured molecular gas mass is 
equivalent to that of the typical SMG which has only a 
weak active nucleus (Alexander et al. 2005). Therefore, 
our observations suggest that the 3C318 system is un- 
dergoing a starburst comparable to that of the typical 
SMG. 

3. CONCLUSIONS 

We have carried out a sensitive search for CO (2-1) 
emission in the z = 1.574 quasar 3C318. We find a 
molecular gas mass of M Ha = (3.0 ± 0.6) x 10 10 M Q . 
This value is typical for SMGs, leading to the conclusion 
that the predominant power source for the exceptionally 
high cool dust luminosity in this quasar is due to star 
formation. 

The implications of these observations go beyond this 
one quasar and provide extra evidence towards the syn- 
chronisation of star formation and jet formation events 
in some radio-loud AGN. Willott et al. (2002) found 
that small (young) radio sources, such as 3C318, are 
more likely to be strong sub-mm sources than large (old) 
radio sources. Additionally, several of the most sub-mm- 
luminous quasars in Willott et al. (2002) have reddened 
optical spectra (like 3C318), indicating significant dust 



in their host galaxies. There is also a tendency for the 
most sub-mm-luminous quasars to have strong Civ as- 
sociated absorption. Baker et al. (2002) presented a 
negative correlation between C iv associated absorption 
equivalent width and radio source size (age). Morganti 
et al. (2005) found that small radio sources have fast 
(~ 1000 km s -1 ) outflows carrying a considerable mass 
of neutral gas. Putting all these lines of evidence to- 
gether points towards a dusty, metal-enriched interstel- 
lar medium with a high rate of star formation in quasars 
with recently triggered radio jets. The larger-scale mas- 
sive gas outflows which follow may be effective at clearing 
out much of this material and suppressing star formation. 
As shown by Morganti et al., such outflows could be en- 
ergetic enough to not only inhibit star formation, but ex- 
pel the remaining gas from the galaxy. Therefore, radio 
sources may be significant contributors to the feedback 
required to explain the observable properties of elliptical 
galaxies and their black holes. 
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